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Abstract 



In this work we consider the role of muons in various URCA processes relevant 

for neutrino emissions in the core region of neutron stars. The calculations 

are done for /3-stable nuclear matter with and without muons. We find muons 

to appear at densities p = 0.15 fm~ 3 , slightly around the saturation density 

for nuclear matter po = 0.16 fm -3 . The direct URCA processes for nucleons 

are forbidden for densities below p = 0.5 fm -3 , however the modified URCA 

processes with muons (n + N^p + N + p+ V^jp + N + p^ n + N + u^), 

where N is a nucleon, result in neutrino emissivities comparable to those from 

(n + iV — > p + N + e + V e , p + N + e — > n + N + u e ). This opens up for further 

possibilities to explain the rapid cooling of neutrons stars. Superconducting 

protons reduce however these emissivities at densities below 0.4 fm -3 . 
PACS number(s): 97.60.Jd 21.65.+f 74.25.Bt 
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Typeset using REVTgX 



The thermal evolution of a neutron star may provide information about the interiors 
of the star, and in recent years much effort has been devoted in measuring neutron star 
temperatures, especially with the Einstein Observatory and ROSAT, see e.g. Ref. 0. The 
main cooling mechanism in the early life of a neutron star is believed to go through neutrino 
emissions in the core of the neutron star The most powerful energy losses are expected 
to be given by the so-called direct URCA mechanism 

n— >p + e + F e , p + e — ► n + u e , (1) 

as discussed recently by several authors ||J4]||]. However, in the outer cores of massive 
neutron stars and in the cores of not too massive neutron stars (M < 1.3 — 1.4M ), the 
direct URCA process is allowed at densities where the momentum conservation k F < k F + k F 
is fulfilled[]. This happens only at densities p several times the nuclear matter saturation 
density po = 0.16 fm~ 3 . 

Thus, for long time the dominant processes for neutrino emission have been the so-called 
modified URCA processes first discussed by Chiu and Salpeter || , in which the two reactions 

n + n— >p + n + e + V e , p + n + e-^n + n + u e , (2) 

occur in equal numbers. These reactions are just the usual processes of neutron /3-decay 
and electron capture on protons of Eq. ([I]), with the addition of an extra bystander neutron. 
They produce neutrino-antineutrino pairs, but leave the composition of matter constant 
on average. Eq. (§) is referred to as the neutron branch of the modified URCA process. 
Another branch is the proton branch 

n + p->p + p + e + V e , p + p + e^n + p + i/ e , (3) 

pointed out by Itoh and Tsuneto ]7J and recently reanalyzed by Yakovlev and Levenfish || . 

throughout this work we will reserve letters n, p, e and p to neutrons, protons, electrons and 
muons respectively. Also, we will set h = c = 1. 
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The latter authors showed that this process is as efficient as Eq. (Q). Similarly, at higher 
densities, if muons are present we have the processes 

n — > p + n + V^, p + /i — > n + Ufj,, (4) 

n + n — > n + p + H + Vp, n + p + fi—^n + n + u^, (5) 

and 

n + p — ► p + p + fM + Vp, p + p + fi->n + p+is^, (6) 

in addition one also has the possibility of neutrino-pair bremsstrahlung, processes with 
baryons more massive than the nucleon participating, such as isobars or hyperons ||||, or 
neutrino emission from more exotic states like pion and kaon condensates PlJlOf or quark 
matter [|TT|— pT^] . Actually, Prakash et al. [§J showed that the hyperon direct URCA processes 



gave a considerable contribution to the emissivity, without invoking exotic states or the large 
proton fractions needed in Eq. ([I]). 

The scope of this Letter is to give an estimate of the processes of Eqs. (fj])-© at densi- 
ties corresponding to the outer core of massive neutron stars or the core of not too massive 
neutron stars, performing a self-consistent Brueckner-Hartree-Fock (BHF) calculation for 
/3-stable nuclear matter. We impose the relevant equilibrium conditions, with and with- 
out muons, employing modern meson-exchange potential models for the nucleon-nucleon 
interaction. In addition, we consider also the role of superconducting protons in the core of 
the star, by calculating the pairing gap for protons in the 1 Sq state. The proton supercon- 
ductivity reduces considerably the energy losses in the above reactions ||, and may have 
important consequences for the cooling of young neutron stars. The final outcome of our 
calculations yields then the equation of state, effective masses m* for protons and nucleons, 
their corresponding Fermi momenta kp, the Fermi momenta for electrons and muons and 
the proton pairing gap. These ingredients enter into the calculation of the various reaction 
rates. 



In more detail, our calculational procedure is as follows: 
First we solve self-consistently the BHF equations for the single-particle energies, using a 
G-matrix defined through the Bethe-Brueckner-Goldstone equation as 

G = V + V Q TT G, (7) 

where V is the nucleon-nucleon potential, Q is the Pauli operator which prevents scattering 
into intermediate states prohibited by the Pauli principle, Hq is the unperturbed hamiltonian 
acting on the intermediate states and u is the so-called starting energy, the unperturbed 
energy of the interacting states. Methods to solve this equation are reviewed in Ref. |14| . 



The single-particle energies for state ki (i encompasses all relevant quantum numbers like 
momentum, isospin projection, spin etc.) in nuclear matter are assumed to have the simple 
quadratic form 

where m* is the effective mass. The terms m* and 8, the latter being an effective single- 
particle potential related to the G-matrix, are obtained through the self-consistent BHF 
procedure. The so-called model-space BHF method for the single-particle spectrum has 
been used, see e.g. Refs. fTijl , with a cutoff ku = 3.0 fm -1 . This self-consistency scheme con- 



sists in choosing adequate initial values of the effective mass and 5. The obtained G-matrix 
is in turn used to obtain new values for m* and S for protons and neutrons. This procedure 
continues until these parameters vary little. The BHF equations are solved for different pro- 



ton fractions, using the formalism of Refs. [14| , |15| . The conditions of (3 equilibrium require 
that 

= Hp + He, (9) 

where \i{ is the chemical potential and that charge is conserved 

Pp = Pe, (10) 
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where pi is the particle number density in fm 3 for particle i. If muons are present, we need 
also to satisfy charge conservation 

Pp = Pe + P^, (11) 

and energy conservation through 

He = /V (12) 
The nucleon-nucleon potential is defined by the parameters of the meson-exchange potential 



model of the Bonn group, version A in Table A. 2 of Ref. |^q| . From the derived effective 
interaction one can in turn calculate the energy per nucleon as function of the total density 
14|. This is plotted in Fig. [I] with and without muons. As can be seen from this figure, 



muons yield contributions to the energy per particle (and in Fig. |2| to the pairing gap) at 
densities around 0.15 fm -3 . The presence of muons makes the energy per particle and in 
turn the equation of state softer, yielding smaller masses and larger radii (M = 1.45M and 
R fh 9 km here) for neutron stars. Similar qualitative results are reached if one performs 
a Dirac-Brueckner-Hartree-Fock calculation, as in Ref. [|I7]], although the equation of state 
in that case is stiffer, resulting in larger masses and smaller radii for the star. More details 
and applications to neutron stars will be reported by us in a future work fl8[| . 

The next step is to evaluate the gap equation following the scheme proposed by Anderson 



and Morel |19j and applied to nuclear physics by Baldo et al. |20|. These authors introduced 
an effective interaction Vfc.fe'- This effective interaction sums up all two-particle excitations 
above the cutoff ku- It is defined according to 

Vk,k' = Vk,k> — X! Vk,k" nj7l Vk»,k', (13) 

k">k M ltj k» 

where the quasiparticle energy is given by 



E k = ^(e k -e F y + Al (14) 

Sf being the single-particle energy at the Fermi surface, Vk,k' is the free nucleon-nucleon 
potential in momentum space and A*, is the pairing gap 



A fc = - £ v^^r- (15) 
k'<k M ltjk ' 

For notational economy, we have dropped the subscript % on the single-particle energies. 

In summary, first we obtain the self-consistent BHF single-particle spectrum Ek for 
protons and neutrons in /^-stable matter. This procedure gives the relevant proton and 
neutron fractions at a given total density. Thereafter we solve self-consistently Eqs. ([13]) 
and ( |i5| ) in order to obtain the pairing gap A for protons or neutrons. For states above Jzm, 
the quasiparticle energy of ([H]) is approximated by Ek = {ek — Sf), an approximation found 
to yield satisfactory results in neutron matter pif . Our prescription gives then the pairing 
gap at T = 0. The results for the proton pairing gap with and without muons are displayed 
in Fig. [| as functions of the total baryonic density. 

The reader should notice that our calculation of the pairing gap does not include com- 
plicated many-body contributions like polarization terms. These may further reduce the 
pairing gap and bring in an uncertainty regarding the size, which strongly influences the 
modified URCA processes. Although we have employed the Bonn A potential in Table A. 2 



of Ref. fL6] , the pairing gap obtained with other potentials like the Paris or Bonn B and 
Bonn C is rather similar, results. This is expected since the x Sq channel is determined by the 
central force component only of the nucleon-nucleon interaction, and all modern potentials 
reproduce the phase shifts for this channel. 

With the pairing gap, effective masses and the Fermi momenta for the various particles, 
we first look at the neutrino energy production rate for Eq. (0). The expressions for these 
processes were derived by Friman and Maxwell p2j and read || 



^-8-5xl0-(^) 3 (5)fe) 1/3 T 9 ^, (16) 

in units of ergs cm _3 s _1 where T 9 is the temperature in units of 10 9 K, and according to 
Friman and Maxwell a n describes the momentum transfer dependence of the squared matrix 
element in the Born approximation for the production rate in the neutron branch. Similarly, 
f3 n includes the non-Born corrections and corrections due to the nucleon-nucleon interaction 



not described by one-pion exchange. Friman and Maxwell used a n ~ 1.13 at nuclear 



matter saturation density and (3 n = 0.68. In the results presented below, we will not include 
a n and [3 n . For the reaction of Eq. (|5|), one has to replace p e with p M . For the proton branch 
of Eq. (|3|) we have the approximate equation [g] 

1/3 



g p ^8.5xio 21 (&) (K) M 1 

^ p \nipj \m n J \poJ 



x T$a p [3 p ( 1 - §r ) 9, (17) 



where G = 1 if k F < 3k F + k e F and zero elsewhere. Yakovlev and Levenfish put a n = a p and 
P P = Pn- We will, due to the uncertainty in the determination of these coefficients, omit 
them in our calculations of the reaction rates. For Eq. (|6|) one replaces p e with and k e F 
with k F . 

The reaction rates for the URCA processes are reduced due to the superconducting 
protons. Here we adopt the results from Yakovlev and Levenfish ||, their Eqs. (31) and 
(32) for the neutron branch of Eq. (|l^) and Eqs. (35) and (37) for the proton branch of Eq. 
(|17|). In this work we study proton singlet-superconductivity only, employing the standard 
approximation 



k B T c 



0.5669, (18) 



A(0) 

where Tq is the critical temperature and A(0) is the pairing gap at zero temperature derived 
in our calculations. The critical temperature is then used to obtain the temperature depen- 
dence of the corrections to the neutrino reaction rates due to superconducting protons. To 
achieve that we employ Eq. (23) of ||. 

In this Letter, we present results for the neutrino energy rates at a density p = 0.3 
fm -3 . The critical temperature is Tq = 3.993 x 10 9 K, whereas without muons we have 
Tq = 4.375 x 10 9 . The implications for the final neutrino rates are shown in Fig. ^|, where 
we show the results for the full case with both muons and electrons for the processes of Eqs. 
(0), (^), (|5|) and (^). In addition, we also display the results when there is no reduction 
due to superconducting protons. At the density considered, p = 0.3 fm~ 3 , we see that the 
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processes of Eqs. (||) and are comparable in size to those of Eqs. (§) and @ (on the 
log-scale there is basically no difference). The direct URCA process of Eq. (|T|) is not allowed 
in our /3-stable matter for densities p < 0.5 fm~ 3 , due to momentum conservation. Similarly, 
the direct URCA process with muons is also not allowed at densities p < 0.6 fm~ 3 . Thus, 
Eqs. ([5p and (^) give an additional and important contribution to the neutrino production 
in the core of a neutron star. However, the proton pairing gap is still sizable at densities up 
to 0.4 fm~ 3 , and yields a significant suppression of the modified URCA processes discussed 
here, as seen in Fig. [|. We have omitted any discussion on neutron pairing in the 3 P2 state. 



For this channel we find the pairing gap to be rather small f21fl , less than 0.1 MeV and close 
to that obtained in Ref . |24| . We expect therefore that the major reductions of the neutrino 
rates in the core come from superconducting protons in the 1 Sq state. Moreover, the recent 



reinvestigation of neutrino-pair bremsstrahlung in the crust by Pethick and Thorsson p5 
indicates that this process is much less important for the thermal evolution of neutrons stars 
than suggested by earlier calculations. Thus, combined with the present reduction due to 
superconducting protons, our results may indicate that there is need for other processes 
than those studied here to explain the rapid cooling of neutron stars. The contribution 
from neutrino-pair bremsstrahlung in nucleon-nucleon collisions in the core is also, for most 
temperature ranges relevant for neutron stars, smaller than the contribution from modified 
URCA processes [§]. Possible candidates are then direct URCA processes due to hyperons 
and isobars, as suggested in Ref. 0, or neutrino production through exotic states of matter, 
like kaon or pion condensation [d,|K]] or quark matter [|IT| -|i"3|| . Though, the analysis of Page 



|| indicates that it is hard to discriminate between fast or slow cooling scenarios, though 
in both greement with the observed temperature of Geminga is obtained if baryon 

pairing is present in most, if not all of the core of the star. 

In summary, we have performed a self-consistent calculation for /3-stable neutron matter, 
with and without muons. We have shown that the modified URCA processes (n + TV — > 
p + N + p + V ll ,p + N + p — > n + N '+Ufj), where N is a nucleon, result in neutrino emissivities 
comparable to those from (n + N — > p + N + e + Ts e ,p + N + e —> n + N + u e ). These processes 
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should therefore be accounted for in a scenario for neutron star cooling. 

This work has been supported by the Istituto Trentino di Cultura, Italy, the Research 
Council of Norway and the NorFA (Nordic Academy for Advanced Research). 
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FIGURES 



FIG. 1. Energy per baryon as function of the total baryonic density p in /3-stable neutron 
matter with (dashed line) and without (solid line) muons. 



FIG. 2. Proton pairing energy gap as function of total baryonic density p in /3-stable neutron 
matter with (dashed line) and without (solid line) muons. 
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FIG. 3. Temperature dependence of neutrino energy loss rates in a neutron star core at a total 
baryonic density of 0.3 fm~ 3 with muons and electrons. Solid line represents Eq. (2), dotted line 
is Eq. (5), dashed line is Eq. (3) while the dash-dotted line is the processes of Eq. (6). The 
corresponding results with no pairing are also shown. 
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FIG. 3. Temperature dependence of neutrino energy loss rates in a neutron star core at a total 
baryonic density of 0.3 fm -3 with muons and electrons. Solid line represents Eq. (2), dotted line 
is Eq. (5), dashed line is Eq. (3) while the dash-dotted line is the processes of Eq. (6). The 
corresponding results with no pairing are also shown. 
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FIG. 2. Proton pairing energy gap as function of total baryonic density p in /3-stable neutron 
matter with (dashed line) and without (solid line) muons. 
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FIG. 1. Energy per baryon as function of the total baryonic density p in /3-stable neutron 
matter with (dashed line) and without (solid line) muons. 
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